The 130 GeV Fingerprint of Right-handed Neutrino Dark Matter 
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Recently, an interesting indication for a dark matter signal in the form of a narrow line, or 
maybe two lines and/or an internal bremsstrahlung feature, has been found in analyses of public 
data from the Fermi-LAT satellite detector. As recent analyses have also shown that there is little 
sign of extra contributions to continuum photons, it is natural to investigate leptophilic interacting 
massive particle (LIMP) models. We show that a model of radiatively generated neutrino masses 
may have the properties needed to explain the Fermi-LAT structure around 130 GeV. This model 
was proposed some 10 years ago, and predicted a clearly observable 7-ray signal in the Fermi-LAT 
(then GLAST) detector. Here, we update and improve that analysis, and show as an example that 
a right-handed neutrino of mass 135 GeV should give rise to three conspicuous effects: a broad 
internal bremsstrahlung bump with maximum around 120 GeV, a 27 line around 135 GeV, and a 
Z7 line at 119.6 GeV (neglected in the previous work). These features together give a good fit to 
the 130 GeV structure, given the present energy resolution of the Fermi-LAT data. An attractive 
feature of the model is that the particle physics properties are essentially fixed, once the relic density 
and the mass of the right-handed neutrino dark matter particle have been set. Puzzling features of 
the data at present are a slight displacement of the signal from the galactic center, and a needed 
boost factor of order 5 — 15. This presents interesting challenges for numerical simulations including 
both baryons and dark matter on scales of 100 pc, and perhaps a need to go beyond the simplest 
halo models. With upcoming experiments having better energy resolution, or with future Fermi- 
LAT data, the double-peak structure with a definite predicted ratio of the strengths of the two lines 
and the internal bremsstrahlung feature should be seen, if this model is correct. With the planned 
satellite GAMMA-400, a striking fingerprint of this dark matter candidate should then appear. 

PACS numbers: 95.35.+d, 14.60.St, 95.85.Pw, 95.85.Ry, 98.70.Rz 
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I. INTRODUCTION 

For 7-rays coming from annihilations of dark mat- 
ter particles in the halo, the Fermi-LAT instrument [l[ 
has very successfully delivered bounds that have recently 
started to probe into the parameter space of viable mod- 
els, in line with pre-launch expectations @], in particular 
for dwarf spheroidal galaxies [|| and galaxy clusters 3 . 

Interestingly, there have very recently appeared anal- 
yses of public Fermi-LAT data Q where a rather strong 
tentative dark matter signal has been seen, in a line- 
search from the halo in the vicinity of the galactic cen- 
tre @, 0] There is a feature, which has been interpreted 
as either originating from internal bremsstrahlung 0, Q 
around 150 GeV 0], or a 7-ray line at around 130 
GeV Q. This discovery has alrea dy g ained considerable 
attention from theorists (see, e.g., [ill] and [lU, [l3[ for re- 
cent reviews). The line signal has subsequently been in- 
dependently verified by other analyses of the same public 
data [3 Ha]. In particular, in the analysis [l5[ a signal 
of more than 5cr was claimed, with some evidence also 
for a second line, with energy consistent with annihila- 
tion into Z7, which generally also should exist in many 
models [l||. 

There have already appeared quite a number of pro- 
posed models for this feature aound 130 GeV [ll|, al- 



though the previously much studied supersymmetric 
models now seem disfavoured in most of the supersym- 
metric scenarios, due to the non-appearance of a char- 
acteristic 7-ray continuum that should result from frag- 
mentation of quarks, Z- and VF-bosons or gluons [17J. 
(See, however, [131 ] for some interesting surviving parts 
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of the huge supersymmetric parameter space.) 

Let us now turn to another attractive model from the 
particle physics point of view, namely radiative see-saw 
models for neutrino masses (for a recent review contain- 
ing many of the original references, see [H[). 

Some 10 years ago, an interesting model was proposed 
[T9I I . where a right-handed neutrino of mass less than 
about a TeV plays a crucial role in giving mass to the 
otherwise massless standard model neutrinos through a 
high-order loop mechanism. (This is a version of the Zee 
model [2(|. See also [HIHH for a more recent formulation 
of models with similar phenomenology.) 

In this model, neutrino masses appear only at the three 
loop level, achieved by supplementing in the simplest ver- 
sion the Standard Model (SM) fields with two charged 
singlet scalars S\ and S2 and one right handed neutrino 
Nr. (We will here stick to this simplest model, as we ex- 
pect the 7-ray features studied here depend very little on 
the details of more elaborate, and perhaps more realistic 
schemes, such as elaborated upon in [lj| [2l], [23[). Lep- 
ton number is broken explicitly by including a Majorana 
mass term for the right-handed neutrino, and imposing 
a discrete Z2 symmetry under which the SM fields and 



2 



S\ are neutral but 5*2 and N R transform as 
Z 2 : {S 2 ,N R } — ► {-S 2 ,-N R } , 



(1) 



which forbids Dirac masses for the neutrinos. This gives 
the Lagrangian 



£ = fapL^CiT2LpSi + gaNRS^lan 

+M R N£CN R + V(Si, S 2 ) + h.c. 



(2) 



The potential V(Si, S 2 ) contains in particular a (S1S2) 2 
coupling, and a mild hierarchy of masses M R < Mg x < 
Ms 2 oc 100 GeV - 1 TeV is assumed. Furthermore, the 
Yukawa couplings f a p, g a are of order unity, making N R 
stable in view of the discrete symmetry. Left-handed 
neutrino masses are first induced at three loops [H, [l!| . 
For Ms 2 of the order of a TeV, one finds an effective 
dimension- five effective mass scale of A > 10 9 GeV, giv- 
ing remarkably neutrino masses at the 0.1 cV scale with- 
out involving fundamental mass scales larger than a TeV. 

In this leptonically interacting massive particle model 
(a leptophilic WIMP, named LIMP in the paper by Baltz 
and Bergstrom [24[, BB in the following), N R becomes 
stable and is therefore a natural dark matter (DM) candi- 
date. Through S 2 exchange it couples extremely weakly 
today, due to the low DM velocities in the galactic halo. 
The reason for this is its Majorana nature, which means 
that the couplings in the S wave are proportional to the 
mass of charged leptons in the final state (25|. In the 
early universe, on the other hand, P-wave annihilation 
was important and set the relic density QjjMh 2 ~ 0.11 
(the scaled Hubble constant h ~ 0.7). It was noted in 
the original proposal [l9| that this would be a fine DM 
candidate, but thought to be essentially undetectable, as 
the rate of direct detection through scattering on nucle- 
ons would be very small. Similarly, indirect detection 
through neutrinos from the Earth or the Sun will not 
be possible, as the cross section for capture is negligibly 
small. 

However, in BB [24 |. it was pointed out that this can- 
didate, on the contrary, has excellent detection prob- 
ability in indirect detection through the annihilation 
NrN r — S- 27 or through the internal brcmsstrahlung pro- 
cess N R N R —> (the last process with its surprising 
avoidance of helicity suppression, and the relation to the 
27 line was first discovered in Q in the context of super- 
symmetric models for dark matter, with refinements over 

the last few years @, H IM E3 ) ■ Thus ' this model has the 
intriguing property of connecting the dark matter prob- 
lem with that of neutrinos masses. The unique window is 
through 7-ray detection of annihilation of the typical fin- 
gerprint of this model: an internal brcmsstrahlung broad 
enhancement with maximum near 90 % of the N R mass, 
supplemented with two 7-rays lines from the 77 and Z"f 
final states, where the 77 line will appear at E y = m R 
and the -Z7 line at E 1 = m R (l — m% / (4m R )) , as dictated 
by energy and momentum conservation. Both lines are 
intrinsically very narrow, meaning that in practice the 
energy resolution of the detector will be decisive for their 



detection. (Note that for Majorana particles, the anni- 
hilation in the S-wave takes place from a pseudoscalar 
state, and therefore no if 7 line, with if the Higgs boson, 
is expected, as a — > radiative transition is forbidden.) 

In BB, Zj electroweak mixing and therefore the Z7 
annihilation channel was neglected (and in @ it was not 
kinematically allowed). That this model has a very good 
chance of giving an observable structure in Fcrmi-LAT is 
in fact strengthened by a more complete analysis, which 
will be briefly presented here. We will find a striking 
agreement of the shape of the structure with the (ad- 
mittedly still scarce) Fermi-LAT data analyzed in Q. 1 
A very interesting aspect, which contrasts with super- 
symmetric DM, is the essential lack of freedom of choice 
of parameters and therefore the predictive power of the 
model. 

The cross determining the relic N R density can be writ- 
ten SI Si 



av (N R N R ->■ £ + £ 



327rm|j (1 + Pf 



1 + / 4 



2 2 1 
m R v + 



(3) 



with rri£ the charged lepton mass, mj = fm R the S 2 
mass and m R the LIMP (N R ) mass. For / not too much 
larger than unity, the factor in square brackets is typi- 
cally close to 0.5, so that for m R not larger than a TeV, 
the P-wave cross section (the term proportional to v 2 in 
Eq. (J3j> ) determines the freeze-out temperature Tf and 
thus the relic abundance. Typically, Tf/m R ~ 1/20, and 
thus (v 2 ) f = 6Tf/m R ~ 0.3. The requirement that the 
LIMP be the dark matter particle, with relic abundance 
£l R h 2 w 0.11, then fixes the freeze-out cross section 



— e,//,r 



av) f 



E^ 4 (i + / 4 ) 

2907rm2j(l + f 2 f 



3 x 10" 26 cm 3 s" 1 . 



(4) 

For a given N R of mass m R , this sets the normalization 
of the combination + / 4 )/( TO fl(l + / 2 ) 2 ) which 

governs many annihilation processes. The total annihi- 
lation rate into lepton pairs at rest (putting v = 
in Eq. (J3J) ) is for each flavour (assuming flavour universal 
couplings) 



(crv) v „o ~ 10 



-25 



3 -1 

cm s 



(5) 



It was noted [19j, |24j that in order that the couplings 
gi not be larger than 1, for nearly degenerate N R and 



There are many details of the model that are worth studying 
more closely, and that may move the results up or down by a fac- 
tor of two or so (see for example [2lll where a second righthanded 
neutrino is shown to be needed) - a small uncertainty compared 
to many astrophysical unknowns, such as the DM density distri- 
bution near the galactic centre. 
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S2, fnR cannot be too much heavier than 1 TeV, but 
could be significantly lighter. We now thus insert the 
value found in of around 130 GeV for tur. Actually 
with rrifj = 135 GeV we find a slightly better fit, and 
we will use that value in the following. 2 Assuming near 
degeneracy, but choosing for simplicity / = 1.1, to avoid 
complications of possible coannihilations (this could of 
course be relaxed in a more refined treatment), we find a 
value of g ~ 0.52, where we for simplicity assume flavor 
universality, i.e., g P = g^ = g T = g (see [2l| for a more 
realistic ansatz in a more elaborate model) . LIMPs with 
these values of the parameters have a very small annihi- 
lation rate for the e + e~ and /U + /i~ channels. Even for 
t + t~, the helicity suppression is of the order of 10~ 4 . 
At that level, the P-wave term in the galactic halo also 
may contribute a small amount. We neglect that here, 
but include continuum 7-rays from decays of directly pro- 
duced r leptons in our Figures. They turn out to be not 
very important, however. Also, there is a contribution 
to 7-rays from Z decays produced in the similar 'gauges- 
trahlung' process with the final state £ + £~Z which also 
is surprisingly large [2H |27j . (There is no corresponding 
VF-strahlung for these singlet fields.) This contribution 
from Z decays is also included in the numerical estimates, 
but also turn out to be inessential. 

In fact, predictions for a 100 GeV LIMP were made 
in BB (see Fig. 2 (a) in [24|), and shown to be possibly 
quite spectacular. (This is in contrast with the predic- 
tions for positrons from the same process which need a 
"boost factor" of more than several thousand to match 
the anomalous positron ratio of which there were indica- 
tions already at that time (24|.) Thus there is a reason 
for scrutinizing present Fcrmi-LAT data to see whether 
the 135 GeV LIMP fits the spectral feature recently dis- 
covered. 3 Indeed, we will find that this model, especially 
when augmented with the Zj line (which was left out 
in BB) and having both that line, the 27 line and the 
internal brcmsstrahlung contribution fits the Fermi-LAT 
feature surprisingly well. We will compare with the origi- 
nal analysis in [7j and choose for convenience the Einasto 
profile Reg. 3 "SOURCE" events in the parlance of that 
reference. 4 

We note that we are not making a "best fit" analysis 



2 One should be aware that there is an up to a few GeV uncertainty 
of the exact value, depending on the present low statistics and 
also factors like the calibration of the Fcrmi-LAT instrument. 
The relative location of the three 7-ray features is however fixed 
once the mass is set. 

3 This feature often is referred to as the "130 GeV line" . However, 
as we will see, in the present model it is a rather complicated 
composite structure, which is why we prefer to call it the "130 
GeV fingerprint", which for thr = 135 GeV consists of a 135 
GeV 77 line, a 119.6 GeV Z7 line and an internal bremsstrahlung 
continuum, from £+£"7. 

4 This is the preferred density profile, see Fig. 3 in a remarkable 
recent analysis [l3T | - actually, the first to use this tentative signal 
to constrain details of the halo distribution. 



as this is not very meaningful given the presently rather 
scarce data. The point is that with more and better data 
in the future, the 7-ray "fingerprint" of this model should 
present itself, if this model for the dark matter is correct. 

To make sure that continuum 7s are not overproduced 
in Z decays from the Zj final state, we include this 
contribution using the freely available DarkSUSY code 
[28| . but also this flux is quite small. In recent versions 
of the DarkSUSY code all three processes: internal 
bremsstrahlung, 27 and Zj are included and well doc- 
umented in the original literature @, d, [29[ , [3(1 [3l[ and 
[32^ . respectively, so here we only show and discuss the 
results. We note that all three processes have to exist in 
this type of model, and have to be included for consis- 
tency. 5 

Normalizing to the correct relic density, the internal 
bremsstrahlung contributes around 6.7 ■ 10 -29 cm 3 s~ 1 to 
the annihilation of these slow LIMPs in the halo. The 
bulk forms a broad bump around E 7 ~ 0.9m/j (see 
Fig. 1). The 27 line gives 2.3 • 10 -29 cm 3 s _1 (includ- 
ing a factor 2 for the two photons), and the Z7 line has 
a (av) v ^o of 1.0 • 10~ 29 cm 3 s _1 . The combined rate in 
the galactic halo is thus 

(H^o = l-0-10- 28 cm 3 s-\ (6) 

a value not to far from the estimate of the experimentally 
observed effect in of (~ 12.7 ± 5.7l||)- lO-^cmV 1 
(for the favoured Einasto profile). In Q, this was in- 
terpreted only in terms of one 7 line, but the limited 
energy resolution of the Fermi-LAT instrument, around 
10 % FWHM (HI, is only barely enough to resolve the 
structure into its separate components (see Fig.[T]). Note, 
however, that [l5[ favour, although not very strongly, a 
two-line structure, especially when selecting data in the 
detector from directions where the photons pass more 
matter in the calorimeter and the energy resolution is 
therefore improved by a factor around two [Hj]. We will 
soon return to our prediction for an instrument with this 
better energy resolution. 

In Fig. [T]is shown the total 7-ray differential flux and 
the respective contributing process. The results are given 
for the particle physics factor d{N 1 av) / dE^ for slow an- 
nihilations such as in the galactic halo, i.e., in units of 
cm 3 s _1 GeV _1 , and are very robust around this mass 
range, as the normalization is set by the relic density 
constraint. 

In Fig. [5] a comparison with of these results is made 
with the Fermi-LAT data as analyzed by Weniger (tJ, 
multiplied by E^, and and arbitrary normalization con- 
stant (of order 10) to match the overall size of the ef- 



5 It would be important if a new independent calculation of the 
loop-induced Z-y rate could be made, as the results of the two 
existing calculations [30I , I33I agree well in models where interme- 
diate W bosons dominate, but in the "bino-like" case of relevance 
here, there seem to be some discrepancies [33L 
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Gamma-rays from N model N Dark Matter prediction for y flux 




FIG. 1: The differential photon spectrum for the process 
NrNr — > ^ + i"~7, 77 and Z7, smeared with the present Fermi- 
LAT energy resolution, AE/E ~ O.f . The total spectrum is 
given by the black solid line, the internal bremsstrahlung by 
the dashed red line, the smeared 77 line by the green dash- 
dotted line and the smeared Zj line by the blue double-dash- 
dotted line. At the lower left corner the small contributions 
from the t + t~ final state as well as that from Z decays can 
se seen. 



feet. The reason for this "boost factor" is presently un- 
known, but in this model it has to be explained by astro- 
physical effects, such as the detailed distribution of dark 
matter near the galactic centre. (So-called Sommerfeld 
enhancement [H| is not expected in this model. There 
may in principle be fine-tuned mechanisms like S-wave 
pscudoscalar resonances or particles with higher electric 
charge running the loop, but we do not employ such exot- 
ica here.) The required boost may be related to another 
puzzle of the signal, which is a displacement from the 
exact galactic centre by around 200 pc. 6 

As can be seen, once the overall strength has been set, 
the agreement with present data is (perhaps fortuitously) 
intriguing. For the rather low average energy resolution 
(10% FWHM) of Fermi-LAT, the double-peak structure 
is barely visible. Improving, however, the resolution by 
a factor of two, which has been done in [l5| (sacrific- 
ing statistics) , by selecting observation angles which give 
large path-lengths of the electromagnetic shower in the 
detector, the twin peak structure is much clearer. This is 
qualitatively in agreement with our results (see Fig. [3]). 7 



6 In fact, a preprint recently appeared [36| where such a displace- 
ment is shown not to be unnatural in simulations of the combined 
baryon and dark matter system. It remains to be seen whether 
the larger than expected density near the emission region can 
also be explained by similar effects. The off-set could perhaps 
also be explained by the low statistics of the tentative signal [37l | . 

7 It has been suggested in [l5ll that Fermi-LAT may change its 



FIG. 2: Comparison of the total 7-ray differential energy re- 
sults (multiplied by E 2 ) for a 135 GeV right-handed neutrino 
dark matter candidate with the Fermi-LAT public data [f|, as 
analyzed by Weniger Q. A simple power- law fit ~ E~ 2S to 
the continuous background has been made, and is also shown. 



In Fig [3] is also shown what one may expect from the 
next generation of 7-ray space detectors with energy reso- 
lution at the one percent level, such as GAMMA-400 [H] 
and D AMPE (see [4(| and references therein) . Given that 
the type of model described here is the correct explana- 
tion of dark matter, the features of the signal would be 
striking. With such an instrument one could start ana- 
lyzing the dark matter halo density distribution in some 
detail. In fact, the property of the fingerprint of this 
model is, besides the two strong lines, the rather broad 
and slightly asymmetric bremsstrahlung bump. The ab- 
sence of this bump would rule out the model. 

The theoretical reason for the necessity of the internal 
bremstrahlung bump and its relation to the line signal is 
quite interesting. It was shown in Q how these features 
are crucial for reproducing the effective axial anomaly in 
these theories (which lack anomalies at the fundamen- 
tal level). In fact, the strength of the 77 line can al- 
most trivially be computed by using the anomaly result 
(|.F| = 1 in [8;]). Also, compact formulas for the internal 
bremsstrahlung contribution can be found there (recently 
checked independently [27| ) ■ The validity of these formu- 
las is more general than for the specific JVr case discussed 
here. The strength of the Zj line is trickier to compute 
due to the non-zero mz, but we can use DarkSUSY (based 



observational search strategy so as to favour these side-ways en- 
tering events. If this can be done technically, this interesting 
proposal could mean that Fermi-LAT may establish the exis- 
tence of this dark matter fingerprint with high confidence over 
the next couple of years. 



5 



on j32[) for its calculation. 
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FIG. 3: The 7-ray differential energy results (multiplied by 
E 2 ) for a 135 GeV right-handed neutrino dark matter candi- 
date are shown, with the present Fermi-LAT energy resolution 
AE/E = 10 % FWHM (solid black line), with a factor of 2 
improvement (red dashed line) and with a future 7-ray in- 
strument, such as GAMMA-400 [H (dash-dotted blue line) 
with resolution at the one percent level. The extrapolated 
power-law ~ E~ 2 ' 6 of the presently measured continuous 7- 
ray background is also shown. 

Of course, one should also be aware of the possibility 
that the feature in the Fermi-LAT may be spurious. Al- 
though it seems unlikely that it is caused by a statistical 
fluctuation, one cannot yet exclude that its origin is in 
some (unknown) instrumental effect. Hence one should 
wait for the results of the analysis of the Fermi-LAT col- 
laboration itself before taking our comparison with data 
too seriously. 

It will also be interesting to sec if data can be inde- 
pendently reproduced by the HESS-II instrument which 
just has seen first light (see [4l|, where also the poten- 
tial of CTA [H and an instrument like GAMMA-400 
[39l | is analyzed). Also, it remains to be seen if the re- 
cent rather statistically shaky indications of a signal from 
galaxy clusters [42| will survive more data. This will 
most immediately come from the Fermi-LAT instrument, 
which will continue to collect an amazing amount of im- 
pressively clean and interesting data over the next few 
years. 

The type of model discussed here has been rather ex- 
tensively studied in the neutrino and particle physics 
communities in recent years, thanks to their ability to 
explain the observed neutrino masses without invoking 
a super-high mass scale. If couplings are not flavour di- 
agonal, dangerous non-observed decays like \i — > ej may 
occur. However, this can be rather easily avoided, with- 
out affecting appreciably the results presented here, by 
invoking discrete symmetries, as explained in [l8l ]. In 



that reference, it is also pointed out that, remarkably, 
in the case of having m,R and mj of the order of 100 
GeV and g ~ 0.5, i.e., precisely the range of values that 
are forced upon us if we want to explain the Fermi-LAT 
structure, also the enigmatic discrepancy of the measured 
(g — 2) M value [43[ with the theoretical prediction [44| has 
a chance to be explained. (See, however, Q for a prob- 
lem in a related generic model.) This is an interesting 
topic for further studies. 

Of course, even if the present indications from the pub- 
lic Fermi-LAT would disappear, this model still gives an 
interesting signal to search for using present and new 
7-ray detectors. The signal would naturally appear be- 
tween 100 GeV and a TeV (perhaps, then, with its canon- 
ically predicted rate, a factor of 10 lower than that indi- 
cated in This energy region will be closely watched 



in the near future 41 1 



One should note that there are interesting very recent 
indications that the dark matter density, at least locally, 
may be a factor of around 3 higher than the generally 
adopted value [4o| . It would be interesting to analyze 
the predicted rates in such models with either a dark 
disk and/or an oblate dark matter halo. This could in 
principle explain the boost of order 10 needed to fit the 
Nr dark matter model to the Fermi-LAT data as ana- 
lyzed by Weniger Q. 

Finally, a note on other ways to test the hypothesis of 
LIMP dark matter. Unfortunately, the scattering cross 
section on nuclei will be very small. This was recently 
shown [46| for an effective description of the 27 vertex. 
Also the radiative scattering process with one photon ex- 
change responsible for scattering on nuclei has a depress- 
ingly low rate (see [I?) for a detailed explanation of the 
reason for this suppression). 

The low scattering rate also means a very low capture 
rate, which hurts detection by neutrinos from the Sun or 
Earth. 

Whether there will be traces of this new sector, only 
interacting leptonically with standard model particles, at 
the LHC remains to be seen. In principle, there could be 
mediators that for instance also would explain the some- 
what too high 27 decay rate of the Higgs boson candidate 
indicated by recent LHC data [H[. It would then, simi- 
larly to in our case, be important to detect and measure 
the rate of H — > Z7 ( |49l - [51| ). which in the standard 
model is closely related to H — > 2j. 



Another possibility for a future linear 



collider 



may be contemplated. The Si and S2 states have unit 
electric charge and would appear as (rather slowly ris- 
ing) enhancements in the total e + e~ cross section, both 
starting around 150 GeV energy per beam. 

There is actually another, challenging, idea to directly 
probe the dark matter nature of the right handed neu- 
trino. The idea is to excite with a strong electron beam 
ambient Nr LIMPs to the nearly degenerate S2, which 
would then immediately decay isotropically back to an 
electron and Nr. This was first suggested and computed 
in the supersymmctric scenario in |47| . and recently re- 
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discovered [52J. It seems that substantial experimental 
development would be needed, however. 

It would also be interesting to work out the possible 7- 
ray signals that could appear from very powerful electron 
jets emanating from centres of active galaxies, along the 
ideas of Ref. [53||, for a 135 GeV Nr. 

Multiwavelength studies of the region around the 
galactic centre, in particular radio data, could be impor- 
tant for these models with positron and electron emission 
(see, e.g.,[54[). It seems, judging from a recent detailed 
study, however, that for the Einasto profile current limits 
are not constraining [55| , but maybe with future experi- 
ments this could give an important cross-check. 

To conclude, for indirect detection of dark matter, the 
excellent Fermi-LAT data - made public to the scien- 
tific community with many excellent analysis tools - has 
opened a very exciting opportunity for studies of detailed 
predictions from models such as the one discussed here. 



Within a few years, the aquired data (and here the ver- 
dict of the Fermi-LAT collaboration itself will be impor- 
tant) could make us accept the model with confidence 
- or disprove it. These are interesting times for dark 
matter studies, indeed. 
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